Formation of discontinuities in flux-saturated degenerate parabolic equations

We endow the nonlinear degenerate parabolic equation used to describe propagation of thermal waves in plasma or in a porous medium, with a mechanism for flux saturation intended to correct the nonphysical gradientflux relations at high gradients. We study both analytically and numerically the resulting equation: ut = [uQ(g(u)x)]x, n > 0, where Q is a bounded increasing function. This model reveals that for n > 1 the motion of the front is controlled by the saturation mechanism and instead of the typical infinite gradients resulting from the linear flux-gradients relations, Q ∼ ux , we obtain a sharp, shock-like front, typically associated with nonlinear hyperbolic phenomena. We prove that if the initial support is compact, independently of the smoothness of the initial datum inside the support, a sharp front discontinuity forms in a finite time, and until then the front does not expand. Mathematics Subject Classification: 35K65, 35B65, 35B40